Competing orders in the Dirac-like electronic structure and the non-hnear sigma 

model with the topological term 
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The Dirac-like electronic structure can host a large number of competing orders in the form of 
mass terms. In particular, two different order parameters can be said to be dual to each other, when 
a static defect in one of them traps a quantum number (or "charge") of the other. We discuss that 
such complementary nature of the pair of the order parameters shows up in their correlation functions 
and dynamical properties when a quantum phase transition is driven by fluctuations of the one of the 
order parameters. Approaching the transition from the disordered (paramagnetic) side, the order 
parameter correlation function at the critical point is reduced, while such fluctuations enhance the 
correlation of the dual order parameter. Such complementary behaviors in the correlation function 
can be used to diagnose the nature of quantum fluctuations that is the driving force of the quantum 
phase transition. 



I. INTRODUCTION 

Electrons in grapheneBi, one atom-thick flake of car- 
bon, effectively behave as relativistic particles, governed 
by the (2-|-l)-diniensional Dirac Hamiltonian: the hop- 
ping of TT orbitals on the honeycomb network of carbon 
atoms results in two Dirac points at the corners of the 
hexagonal Brillouin zone {K±). The Dirac electronic 
structure by now is not an exclusive feature of graphene, 
but also appears on thepSurface of the (3,-^.1 )-dimensictfjal 
Z2 topological insulatorS such as Bi2Se3ul'El and Bi2Te3BLl, 
and in the quasi-ptwo-dimensional organic conductor a- 
(BEDT-TTF)2l3.ffl 

In the Dirac electronic structure, a mass gap can be 
induced in the fermionic energy spectrum by a non- 
vanishing and uniform order parameter. In graphene, 
a variety of order parameters take on a form of such a 
mass term, including, e.g., the antifcrromagnetism, va- 
lence bond solid (VBS), charge density wave (CDW), and 
superconductivity. Some of these mass terms arc of topo- 
logical nature, in the sense that when added they realize 
a topological insulator, such as the quantum Hall effect 
(QHE) or quantum spin Hall effect (QSHE), or a topo- 
logical superconductor. (See, for example, Refs. 10 -18 
and references in Ref. p^ . 

Order parameters (mass terms) proximate to the Dirac 
electronic structure can be classified according to their 
algebraic properties. While they all anticommute with 
the Dirac kinetic term, as they induce a mass gap, they 
can either mutually commute or anticommute with each 
other (See below for more details). A set of mass terms 
which anticommute with each other do not compete, in 
the sense that one can adiabatically changes one order 
into another without closing a gap in the fermionic spec- 
trum. Such mass terms can then be unified into a single 
multi-component order parameter, and each mass can 
be viewed as representing a different direction in the or- 
der parameter space. For example, the three antiferro- 
magnetic orders (Neel orders) and the two different pat- 



terns of the VBS do not compete and, at the level of 
the fermionic spectrum, can be integrated into a five- 
component (0(5)) vector. 

When a topological defect is introduced in such a 
multi-component order parameter, midgap states appear 
in the massive Dirac energy spectrum. For example, for 
the case of the Neel- VBS 5-tuplet, one can create a vor- 
tex, say, in the x and y components of the Neel vector, 
or, in the two components of the VBS order parame- 
ters. A closer look at these topological defects and the 
midgap states reveals an interesting "duality" relation: A 
topological defect in the VBS order parameter supports 
two midgap states which carry a spin quantum number 
(i.e., form a doublet of S = 1/2). On the other hand, 
a topological defect in the Neel vector, again, supports 
two midgap states which can be viewed as forming a dou- 
blet of the "valley pseudospin" degree of freedom. Here, 
by "valley pseudospin" , we mean the degree of freedom 
associated with the two Dirac cones (valleys) in the hon- 
eycomb lattice, and regard two valleys at K± as "pseudo 
spin up" and "pseudo spin down" , respectively. Such du- 
ality relation among pcder parameters can be found also 
in (3+1) dimensions.E^ 

So far, the order parameters have been treated as a 
static background; they are either treated at a mean field 
level, or induced externally. Accordingly, topological de- 
fects are also introduced statically and by an external 
mean. In this paper, we turn our attention to the effects 
of electron correlations and on the dynamics of such order 
parameters. 

While the vanishing density of states in the Dirac spec- 
trum at the half-filling prevents order parameters from 
developing long-range order, when electron-electron in- 
teractions, say, arc sufficiently strong, such orders can 
still be induced. 

Here, an important implication of the duality relation 
between order parameters, found solely at the level of 
the fermionic energy spectrum, is the possibility of "de- 
confined" quantum criticality: I.e., a possible continu- 
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ous (second order) quantum phase transition between 
symmetry unrelated quantum phases, which is forbid^ 
den in the conventional Landau-Ginzburg framework.cll 
The duality relation between different orders, i.e., a topo- 
logical defect in one of the two carry a quantum num- 
ber ("charge") of the other, is a prerequisite for decon- 
fined quantum criticality; with the duality, collapsing one 
phase by condensing defects automatically leads to the 
emergence of the other phase. I.e., they are connected 
through a continuous transition, and do not coexist. 

The deconfincd quantum criticality scenario has been 
investigated so far in the square-lattice quantum antifer- 
romagnet;. Indeed, a recent numerical study of the J-Q 
modc£3n3 observed a vortex-like texture of the VBS or- 
der parameter induced around a deficit of spin in the 
both columnar and plaquette VBS phases. 

In this paper, including sufficiently strong short-range 
electron-electron interactions, we induce a mass gap in 
the Dirac spectrum dynamically (spontaneously), and 
study the structure of the resulting phase diagram and 
critical properties, with eyes on to the scenario of decon- 
fincd quantum criticality. In Sec. ^ we introduce the 
Dirac Hamiltonian together with mass terms that repre- 
sent order parameters of various kinds. The four-fermion 
interactions which result in a spontaneous generation of 
these order parameters are also introduced. One of our 
main results is the anomalous dimensions of the masses 
at the non-trivial critical point, Eqs. ( |l7|) and (|l9| ) in 
Sec. pill The detailed derivation of these a re p resented 
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in Subsecs. Ill B HID 



We conclude in Sec. ^ with our 
speculation on the two-parameter renormalization group 
(RG) flow in the non-linear sigma model (NLcrM) with a 
topological term (0-term) in (2+1) dimensions. 



A. fermion mass terms in graphene 

We now introduce mass terms in the Dirac Hamil- 
tonian ?^o(k). To discuss all possible duality relations 
among order parameters discussed in Ref. ^ in a unified 
fashion, let us start from a set of seven 2^ x 2'^ anticom- 
muting hermitian matrices. 



with ^7 
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(2) 



These matrices form a spinor representation of S0(7). 
Two out of these seven matrices, ^1^2, can be used to form 
a Dirac kinetic Hamiltonian, whereas the remaining five 
matrices can be used as a mass matrix representing an 
order parameter: 



^ — ^i—l ^ — ^1—3 

H = j d^fcl'^H(k)*. 



(3) 



Here TOi=3....j g R represents a five-component order 
parameter; 4* is a fcrmionic field operator which includes 
sublattice, valley, spin, and particle-hole gradings. For 
the cases of our interest, it jturns out we can always 
take ^ to be eight-component.Ej The imaginary-time La- 
grangian is given by 



£ = Vfrt i^d^ + ^^^^ fc^^^ + Y^^^^ j q,. (4) 



B. four-fermion interactions 



II. DIRAC HAMILTONIAN AND ORDER 
PARAMETERS 

We start by describing the Dirac kinetic term of 
graphene Hamiltonian, which is represented by the single 
particle Hamiltonian Ho(k) in momentum space as 

Ho(k) = fc^(T3 (g) Tl ~f fcyCTs T2, (1) 

where 171,2.3 and ti_2.3 are two independent sets of 2 x 2 
Pauli matrices, and we have set the Fermi velocity to be 
unity, vp = 1. In graphene the Pauli matrices T^^y^z act 
on the valley index whereas the Pauli matrices cr^.j/.z act 
on the sublattice index. 

To describe several order parameters, including e.g., 
magnetic and superconducting order parameters, we in- 
troduce two more gradings, one for spin (represented by 
the Pauli matrices 51,2,3), and the other for particle-hole 
(represented by the Pauli matrices /ii,2,3). (For more 
detailed description of the order parameters and the cor- 
responding mass matrices entering in the Dirac Hamilto- 
nian, see Ref. p^ .) 



So far, we have been discussing order parameters rrii 
and their mass matrices (i = 3, . . . , 7) entering in the 
single-particle Dirac Hamiltonian, without asking their 
microscopic origin, i.e. microscopic interactions which 
can generate these orders dynamically. An interacting 
Hamiltonian which gives rise to these orders is: 

-Jd'rj2'-^{^^^a^)\ (5) 

a— 3 

where ga > 0. Such model can be designed on the lat- 
tice with extended Hubbard type interactions. With the 
Hubbard-Stratonovich (HS) transformation, these quar- 
tic interactions can be decoupled into channels, such 
as antifcrromagnctic, superconducting, and VBS orders 
that we have discussed. If the coupling constants ga > 
are large enough, we find a saddle point where these order 
parameters are non-zero. 

While all coupling constants ga can take, in principle, 
different values, below, let us first consider an interacting 
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Hamiltonian with = 5 (a = 3, . . . , 7), 
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2 ^ 



(6) 



This is an 0(5) symmetric analogue of the (thr^ 
dimensional version of the) Gross-Neveu (GN) model.Eil 
The GN model with 0(1) ~ Z2 internal symmetry and 
the chiral GN model with 0(2) ~ U(l) internal symme- 
try have been studied extensively, as a prototype model 
to discuss spontaneous chiral symmetry breaking. The 
corresponding Euclidean Lagrangian is 



£ = f ^ (dr + Cifcx + 6fc,) * - f E (*^^»* 



(7) 



We can successively break this 0(5) symmetry down 
to 0(4), to 0(3), to 0(2) X 0(2), etc. In the following, in 
addition to the 0(5) symmetric model, we consider 0(4) 
and 0(3) symmetric models as well. These models can be 
obtained from Eq. (||) by setting (53,94,5,6,7) = (0,5) (for 
the 0(4) symmetric model) or (53,4,55,6,7) = (0,5) (for 
the 0(3) symmetric model). In these cases, we can treat 
the remaining ^ matrices (i.e. ^3 for the 0(4) symmetric 
model or ^3 and <^4 for the 0(3) symmetric model) as 
mass terms. Schematically, 

beM 

■jd?r Y.i'^Ha^f, (8) 
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where M = 0,{3},{3,4}, for the 0(5)^0(4), and 0(3) 
symmetric models, respectively, and M is the comple- 
ment of M: We have added mass terms X)6ei\f 
which are not generated spontaneously by the interac- 
tions. Rather, they arc entering here as a parameter that 
changes the band structure of the single particle Hamil- 
tonian. We will explore the phase diagram in terms of 
the coupling constant 5 as well as the masses m;, below. 

When the coupUng constant(s) is(are) large enough, 
the system is in an ordered phase. The phase diagram can 
be explored by saddle-point or mean field approximation. 
With the HS transformation. 



If we freeze the dynamics of the HS fields (wa), we get the 
meanfield Hamiltonian discussed before. In the following, 
however, we are interested in the case where the HS fields 
are dynamical. 

The O(A^s) (A^E = 3, 4, 5 is the number of matrices 
introduced below) symmetry of the problem can be seen 



most clearly if we make the following change of basis in 
the fermionic path integral variables: 

70 := -66^5^66, 7i:=-«7oCi, 72 := -j7o6, 
T.a ■■= (a = 1, . . . ,iVs), 

4> := *^7o, tp *, (10) 

wherein the Lagrangian in terms of the new variables is 
given by: 



Y-^^ 2 



(11) 



Here, the summation over space-time index /i = 0, 1, 2 is 
implicit, and we have switched off the masses mj, {b S M) 
momentarily. The merit of this change of variables is that 
it untangles rotations in the order parameter space and 
in the real space: the mass matrices (S^^i jve) and the 
matrices entering in the Dirac kinetic term (7/i=o,...,2) are 
made mutually commuting. 



(12) 



where 7s and Es form S0(3) and SO{N^), respectively. 



H,v = 0, 1,2, 
a,6- l,...,iVs. (13) 



The O(iVs) rotation acting on the fermion fields ip 
can be generated by unitary matrices exp [iX"''] where 
X°-^ = [J^ai^b]- With simultaneously rotating the N-^,- 
component vector Va appropriately, the Lagrangian is left 
invariant. Below, we will denote the dimension of Es and 
7s by Dj: and D^, respectively. In the above construc- 
tion, Ds = Dj = 8, while our calculations presented 
below are valid for any D-^ and Dj. 



C. meanfield phase diagram 

To control the saddle point approximation (the mean 
field theory), and to prepare for the subsequent large- 
Nf expansion, we generalize our Lagrangian to include 
Nf flavors of fermions. In the imaginary time formal- 
ism (the Euclidean signature), the large- A^j generalized 
Lagrangian is then given by 

Nf 

C = 'YtJj,, (j^df, + UaEa + E^ TOfcEfc^ 1p, 

where 7^s are Euclidean gamma matrices. 

We now look for a spatially homogeneous and O(iVs) 
symmetric saddle point solution by setting 



v{x) = \v\n, n ■ n = 1. 



(15) 
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FIG. 1. The mean field phase diagram near the Dirac point. 



When nil, = the self-consistency condition for |f | is 
given by 



DM 



<Pk 



1 



(27r)3 /t2 + |„|2 2Nfg 



= 0. 



(16) 



Here we have introduced the uhra-violet (UV) cutoff A 
by hand. 

The resulting mean-field phase diagram is depicted in 
Fig. |l], wherein we also include to;, (6 S M): For suf- 
ficiently small g <C l/(Af/A), there are two disordered 
phases where fcrmion fields are completely gapped by 
the mass m^gA/ 7^ 0. They are separated by a phase 
boundary mf, = 0. In particular, when mi, represents the 
Kane-Mele mass term (QSHE mass term) it is a quan- 
tum phase transition separating the trivial band insula- 
tor and the topological (QSH) insulator. For sufficiently 
large g ^ l/(iVjA), 0{N^) symmetry is spontaneously 
broken. The arrows in Fig. ^ indicates the infra-red 
(IR) renormalization group (RG) flow. Just at the non- 
interacting Dirac point, four fermion interactions are ir- 
relevant, whereas the mass terms are relevant, from the 
power-counting. The transition between two disordered 
phases is then described by the non-interacting Dirac 
point. The nature of the phase boundary between the 
ordered phase and disordered phases is more difficult to 
discuss. The critical point at mi, = with gNfA 7^ can 
be nevertheless accessible with the large- A^/ or e expan- 
sion (4 — e expansion) . as we explore in the next section. 



III. LARGE-TV/ AND e EXPANSIONS 

A. summary of results 

In this section, wc focus on the nature of the critical 
point separating the ordered phase and the semi-metallic 
Dirac phase located at = with gNfA ^ 0. This IR 
unstable fixed point will be called the GN fixed point. 
While it is not perturbatively accessible in the coupling 
constant g, it can be studied eitheE by the large- A^/ ex- 
pansion or by the e expansion E§cZI The problem of the 
(2-1-1 )-dimcnsional Dirac fermions interacting via short- 
range as well as long-range (Coulombic as well as gauge) 



interactions has a long history; in particular, the re- 
cent fabriefition of graphene, and in particular suspended 
graphcnec3, motivated many papers to revisit this prob- 
lem (see for example, Ref. p9|) . 

One of our main focuses in this paper is the effects of 
four-fcrmion interactions on the duality relation among 
order parameters. Of particular interest are the scal- 
ing dimensions of the mass bilinears at the GN fixed 
point associated to the order parameters; we can com- 
pare the scaling dimensions of the mass ipip and the mass 
ip'^'ip, where the mass matrix S corresponds to one of 
the masses described by (S € 7°^b with h G M), 
and S ^ {^a\a=i,...,Ns- 1-6., S is an identity matrix 
in the Dirac indices, and anticommutes with Sq for all 
a ~ 1,...,Ny: (i.e., the mass ip'E'il) is dual to 
whereas the mass V'V' commutes with -0 

To leading order in the large- A^/ expansion, the anoma- 
lous dimension of the fermion field tp (j]il>)7 the mass bi- 
linear 'lp^p (jliiip)^ ^i^d the mass bilinear -0 !]■(/; {r]^-^^), are 
given by 





8A^s 


'70 = 






16 Ns 




Sn^N 




16 N^ 




3n^N 



64 TVs 
_32iVs 



where we have introduced 



A^ := D.,D^Nf. 



(17) 



(18) 



Alternatively, we can also lo a k .at the anomalous dimen- 
sions using the e expansionEJl^ (i-e., in d = 4 — e dimen- 
sions where physical dimension corresponds to e ~ 1). 
To leading order in e. 



1 



eA^y 



2N + 4-N^ 



A^ + 4 - A^E 

eNj: 
N + i-N^ 



(1 + 2) = 
(1-2) = 



3eA^y 



A^ + 4 - A^s ■ 
-eN^ 

TV + 4 - A^s ' 



(19) 



The calculational details of Eqs. ( |l9|) and (^) can be 
found in the next subsection. 

While the results in the large-A^/ and the e expan- 
sions do not match numerically, the signs of the anoma- 
lous dimensions agree in both calculations. In particu- 
lar observe that 77^^ and rj^-^^ have opposite sign; the 
positive/negative sign in these anomalous dimensions di- 
rectly comes from anticommutation/commutation rela- 
tions of the masses and order parameters. I.e., the pos- 
itive/negative sign is a direct manifestation of the du- 
ality. To be more precise, observed that in 77^^ and 
Ty^Si/" s-bove, two contributions are separately displayed 
(as seen in "(1 ± 3)" and "(1 ± 2)" in Eqs. (|l|) and 
(p7|), respectively). The first terms in and ij^^^ 



5 



-o +-0 o 

FIG. 2. The boson self-energy in the large-TV/ expansion 
obtained by summing fermion bubbles, where full lines repre- 
sent the fermion propagator, and broken lines represent the 
boson propagator. 



come from the renormalization of the single fermion field 
(= 2r]^) and do not know such duality relation, while 
the second terms have opposite sign ("±3" or "±2") for 
Tj^^ and rj^-^^ which originates from the anticommuta- 
tion/commutation relations of the masses and order pa- 
rameters. Such anticommutation/commutation relations 
also control presence/absence of topological terms when 
one integrates the fermionic degrees of freedom to derive 
the non- linear sigma model (NLctM) (see below). 




FIG. 3. The fermion self-energy diagram where the broken 
line represents the boson propagator, Eq. (pi|). 




FIG. 4. The vertex correction diagram relevant for the 
renormalization of the mass bilinear iptp or ijiTii}), where the 
square represents either the insertion of the operator vipip or 



B. Large- TV/ expansion 

We now present some calculational details of the large- 
Nf and e-expansions. In Sees. [II B and [II C, we first col- 
lect diagrammatic calculations for renormalization con- 
stants, in the large- iV/ and e-expans ions, respectively. 
These data will be used in Subsec. [II D to compute 



anomalous dimensions of the fermion field and mass bi- 
linears. 

The bare Euclidean action of the 0{N^) GN model in 
terms of the fermionic [ij)) and bosonic fields {vg) is Eq. 



N, 



1=1 



^ a=l 



(20) 



where we work in d = 3 space-time dimensions, ^ = 

a. Boson self-energy To leading order in the large- 
Nf expansion, the boson {va field) self-energy is obtained 
by summing over fermionic bubble diagrams (Fig. |^), 
leading to the boson propagator 



«(k)^'6(k')) 



16 5ab 5] 



k.k' 



N 



(21) 



sponding integral is: 

(^(q)V'(q))-' 



i4 + 



(27r)3p 7V|q-k| 
A' 

7\ 



(22) 



where ^ ^ A:^7^, and S^J^ = 16iVs/(6^2^). In Eq. (|2|), 
and from now on, we suppress the flavor index, since cor- 
relation functions that consists of several fermionic oper- 
ators are non-zero only when all fermion fields carry the 
same fiavor index: E.g., ('0(q)i/'(q))"^ = (V'<,(q)V't(q))~^ 
in Eq. (p2| ). 

c. Tptp bilinear vertex correction To compute the 
anomalous dimension of the fermion bilinear ipip^ we add 
a mass bilinear term vipil) to the action, where v is the 
corresponding "mass". We also introduce a multiplica- 
five renormalization constant for ly, Ai, , in the large- 
Nf expansion. The vertex diagram corresponding to the 
renormalization of the mass v is Fig. ^ 

A^.log(| 

d^q i{i + l^) IQ^gSab^b i{i - ^) 
(2^)3(g + fc)2 Nq [q-kf 



at the non-trivial critical point g ~ g^. 

To compute the anomalous dimensions of the fermionic 
field, as well as the anomalous dimensions of mass terms 
of different kinds, at the non-trivial critical point, there 
are three diagrams to be considered: Figs. 

b. Fermion self-energy The diagram corresponding 
to the Fermion self-energy is shown in Fig. ^. The corre- 



d. tpTitp bilinear vertex correction We now add a 
term of the form pip'Eip to the action. Here, the mass 
matrix E corresponds to one of the masses described by 
rub (i.e. S € -f°^b with b e M), and E ^ {Ea}a=i,...,Wj:: 
E is an identity matrix in the Dirac indices, and anti- 
commutes with E^ for all a = 1, . . . ,N-^. The diagram 
that determines the renormalization of pipTitp is similar 
to the one given above (Fig. H) except for the presence of 
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FIG. 5. The vertex diagram for A coupling. 



S matrix: 



N 



Ap V log 



A 



= P 



(2^)3(q + fc)2 Nq {q-kf 



C e expansion 

The bare Euclidean action S = J (fir C of the Gross- 
Neveu- Yukawa model, in terms of the fermionic [tp) and 
bosonic fields (wa), is given by 



a=l 



.a=l 



(25) 



We now work in d = 4 — e space-time dimensions. The 
bare propagator for the bosonic field Va is 



«(k)z;,(k')) = 



5ab5\ 



k,k' 



(26) 



To one-loop, we need to consider, for the part of the 
theory that involves only the bosonic field (wa), the di- 
agrams that appear in renormalizing the 0^ theory. In 
addition, the vertex correction for coupling A, and the 
boson self energy correction by a fermion loop need to be 
calculated. 

e. Vertex correction for coupling A The vertex cor- 
rection of the coupling of the fermionic and bosonic fields 
A is given by the following diagram: 



e 
cf^q 



{27ry [(q + p^)2 + rn^] q^q - ky 



Sa(-iVs + 2)i. 



(47r) 



(27) 



/. Boson self energy The one-loop boson self energy 
including a fermion loop is given by 



«(p)^;fc(p))-l 



cfiq 



SabiV^ 



(2-kY 
2 



tr 



d'^q 



:q + p)' 

5ab 



{2tiY q2((?+p)2 



N p2 

(47r)2 T 



5ab 



(28) 



g. Fermion self-energy The diagram corresponding 
to the Fermion self-energy has the same structure as in 
the large Nf expansion (Fig. ||) and the corresponding 
integral is given by 

(^(k)^(k))-i 



iJ/i + A 



2 / ii ^aSab^'' 



{2n)d q^ (k - q)2 



1 + , 



(29) 



where = X^Ns/i'inY. 

From these three diagrams, we read off the renormal- 
ization conditions for the wavefunction renormalizations 
and for the fermionic (-0) and bosonic (va) fields, 
respectively, and for the coupling A as 



•2 

Zli^ + A - — 7^^^ ~ finite, 



(4^) 
,1 N 



e (47r) 



finite, 



2A3 1 



(30) 
(31) 



m-Ao(Z:)^Z^-^-(-iVs + 2) = A, (32) 

where we have introduced an arbitrary mass scale m oc 
-y/r, and Aq is the bare coupling which does not flow with 
the mass scale. With the minimal subtraction schemcEj, 
these conditions can be used to derive, to one-loop, Z^ 
and the beta function for A, 



Z^, = 1 - 



1 N 



e (47r)2' 

(3x2 := -m Aa2 
om 



(33) 



-6A2 + 



2A4 



(47r)2e 

and also the critical coupling 

y2_ i 



[n/2 + 4 - iVs) 



+ 4 - TVj: 



(34) 



Below, we will compute the anomalous dimensions of 
the fermion field and the mass bilinears at the non-trivial 
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Es s 



Villi, 







1 SiVs 


1 e ATe 


6 TX^N 


O A 7 1 A AT 

^ A* + 4 — A* E 


8iVE 


e A^E 




iV + 4 — A^E 




e A^E 




iV + 4 — A^E 


-4 8iVE 


-3 e A^E 


Q aV 


O '\T 1 A AT 

^ Af + 4 — A^E 


2 SA^E 


1 e ATe 




O A T I A AT 

^ A^ + 4 — A^E 


8 8AfE 


3e A^E 


3 7r27V 


iV + 4 - A^E 


-4 8A'e 


e Ns 


3 TT^iV 


N + 4-Nj: 



TABLE I. Anomalous dimensions derived from the large- A'/ 
and e expansions. 



critical coupling (A,ii) = {X*,u*) ^ (0,0). This fixed 
point in the Gross-Ncvcu-Yuka"wa. model is smoothly con- 
nected to the GN fixed point.E3'E3 As in the large-TV/ 
expansion, we will consider the following three diagrams. 

h. tjjip bilinear vertex correction To compute the 
anomalous dimensions, we add the mass term vtpip to 
the action, which leads to the vertex diagram of type 
shown in Fig. ^. It determines, to leading order in e, the 
renormalization constant for the mass term i/'i/jtjj as 



1 



^aSabi^H 



^''""e ""^^ J (2^)rf[(q-fk)2 + m2]g2(g)2 
2A2A^e 1 



(47r)2 



(35) 



i. ip'E'tp bilinear vertex correction Similarly, by 
adding the mass term pipTi^ to the action, the renormal- 
ization constant for the mass term -^pHil) is determined 
from the vertex diagram of type Fig. H: 



^ e 



(47r) 



{2t:Y [(g + p)2 + m2] q2(^)2 
2A2 1 



2A2 1 



(36) 



(47r)2 



D. 



Renormalization conditions and anomalous 
dimensions 



With diagramatics in Sees. Ill B and [II C in hand, we 
can proceed for calculations of anomalous dimensions us- 
ing minimal subtraction schemetZI. The renormalization 



conditions are given by 

ii) {Z^, + Y.-^, Div) = finite, 

Z'^Po + A^ Div p^p, 



(37) 
(38) 
(39) 



where s = Nf or s ~ e for the large- A'^^ and e expansions; 
Div represents the divergence, Div = log(A/fc) for the 
large-Nf expansion and Div = 1/e for the e expansions; 
vq and po are the bare masses. We can readily read off 
the field renormalization Z^, for the fermion field and its 
anomalous dimension 77^,, as 



V4, 



S^Div, 



(40) 



The dimension do of an operator O at a critical point 
is defined from the power-law decay of its correlation 
function as (e'(r)O(O)) oc |r|-2''o. In d = {2 + 1) di- 
mensions, the engineering dimensions of the mass terms 



ipi/j and ipY,iJj are 2, = d^J^^^ = 2. Accordingly, 

the engineering dimensions of the masses ly and p are 



7(0) 



l,d 



(0) 



1 and d%^ ^ d df^ = 



iJi 



From the renormalization conditions (|39D, to leading 
order in Nf or in e. and in d = {2 + 1) dimensions, the 
anomalous part of the dimensions of the masses v and p 
can be read off as 



77^ = A^-S;, 
7?; = A;-S;. 



(41) 
(42) 



The total dimensions of the masses v and p are then given 
by dl ~ Eind d^ ~ l+?7p. In turn, the dimensions of 

the fermion biUnears are then given by d^^ = d — df, and 
^IsV' " '^^^p "^ith the corresponding anomalous dimen- 
sions given by 77^^ = 2c?^^ — 2d and 7 
These results are summarized in Table I 



IV. DISCUSSION 

summary of results Two different order parameters, vi 
and V2 are said to be dual to each other when a static de- 
fect in one of them traps a quantum number (or "charge" ) 
of the other. In this paper, we have discussed the inter- 
action effects on such a pair of order parameters. The 
complementary nature of the pair of the order parame- 
ters shows up also in their dynamical properties (correla- 
tion functions) in the following sense: When a quantum 
phase transition is driven by fluctuations in one of the or- 
der parameter, vi, say, approaching the transition from 
the disordered (paramagnetic) side, the order parameter 
correlation at the critical point is reduced. On the other 
hand, such fluctuations enhance the correlation of the 
order parameter V2 , which is dual to vi . 

If we can experimentally monitor, at a given quan- 
tum critical point, the correlation functions of a pair of 
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such order parameters, the complementary relationship 
(duality) can be revealed by extracting their anomalous 
dimensions, which can, at the same time, be helpful to 
identify the main driving forces of the transition, as dis- 
cussed above. For example, for the Neel and VBS order 
parameters which can compete with each other around 
a quantum critical point, the neutron scattering exper- 
iment can detect the Neel order, whereas the scanning 
tunneling spectroscopy measurement can extract some 
information on the VBS order. Another probe which 
can potentially be used to detect the correlation of the 
VBS order parameter is its coupling to phonons. As a 
consequence of such coupling, information on the VBS 
correlation can be extracted by looking-at the imaginary 
part of the phonon response functior£3 using, say, X- 
ray scattering. Similar method has been used to identify 
fluctuations in the dimer order parametet-in spin-chains 
which go through spin-Peierls transition. 

We, again, emphasize that the type of competitions 
among order parameters discussed in this paper has a 
topological origin, rather than some accidental energetic 
reasons, and hence is largely independent of microscopic 
details of the system. Whether or not such duality rela- 
tionship among correlation functions holds beyond Dirac 
systems that we discussed in terms of the large- A^y and 
e expansions, should further be studied in future. For 
example, numerical study in the J — Q modc]E3 would be 
promising in this regard. 

description in terms of the NLaM To further illustrate 
the topological origin behind the duality, let us seek for 
an alternative description of fluctuating order parame- 
ters. Instead of taking the non-interacting Dirac fermions 
as a starting point, consider the system which is deep in- 
side the ordered phase (the large-gNf region in Fig. ||). 
It is then reasonable to employ a description in terms of 
the fluctuating but non-vanishing order parameter. The 
dynamics of the order parameter may be described by 
the 0{N-s) NLcrM whose kinetic term is given by, the 
following imaginary time action 

SNLaM = J J (fixdfj,nad^na (43) 

where "a = 1 [see Eq. (|l5|)1, and t is the coupling 

constant of the NLcrM. Such NLcrM description can be 
derived, inside the ordered phase, by integrating over 
gapped fermions ip-the presence of slowly varying order 
parameter field vil5. The NLcrM is in the weak coupling 
regime (small t) when g (in Eq. |], say) is large, while it is 
in the strong coupling regime (large t) when g approaches 
the critical coupling (7* from the disordered side. 

In addition to the NLcrM kinetic term, the integra- 
tion over fermions gives rise to either the Wess-Zumino- 
Witten (WZW)-term, 0-term, or the Berry phase term 
in the NLcrM (for alternative derivation of topological 
terms sec Ref. |3^). For example, when we are in the 
0(4) ordered phase, the NLcrM is supplemented with the 



0-term, 

Se = / d XEabcd nadrUbOxncdyUd. (44) 

Here, the 9 aagie can be computed from the microscopic 
Dirac modelEil, 

6* 9 1 

- = 1 - - cos 93 -f- - cos 3(^, (45) 

TT O O 

where ip is determined from the ratio between the ampli- 
tude of the 0(4) order parameter (as defined in Sec. |IIB| ) 
and the mass term, teaup = — \v\/m3: When fermions are 
gapped solely because of the 0(4) order parameter (i.e. 
ma = 0) D = \v\n, which develops a finite expectation 
value in the 0(4) ordered phase, 9 = tt. On the other 
hand, for example, if we add a mass term m^ip'S^'tp in the 
ordered phase, which anticommutes with the 0(4) order 
parameter Vatp'^aip: the 0-term deviates from = tt. 

Folloming the spirit of the Goldstone-Wilczek 
formulaEj, when the order parameters are treated as 
a static background, we can see that it is the ^-term 
which implements the duality within the NLcrM field 
theory. When this operator acquires an expectation 
value, the VBS order results. In this manner the NLcrM 
field theory contains the ingredients describing both 
the ordered and paramagnetic phases. This is in fact 
in line with the known fact that the Berry phase term 
in the (2+1) dimensional 0(3) NLcrM secretly encodes 
the VBS order, and that in the paramagnetic phase-the. 
presence of the Berry phases leads to the VBS order. Elltll 

speculation on the criticality and the RG flow of the 
NLaM with 9 term 

Having established the duality relation within the or- 
dered phase of the NLctM, it is now interesting to ask 
what would happen when we approach the strongly cou- 
pled region where the order parameter fluctuations are 
strong. In particular, how does the RG flow of the 0(4) 
NLcrM with 6'-term look like? The relevance of the topo- 
logical terms (0-term as well as other topological terms 
such as the Wess-Zumino-Witten term) in the NLcrM to 
the deconfined quantum criticality has been discussed in 
Refs. and |39[ (See also a recent a paper |4C| ). 

Unfortunately, the role of the topological term on 
the quantum criticality in (2-|-l)-dimensions is not well- 
understood. This situation should be contrasted with 
the state of our understanding on the 6'-term in the 0(3) 
NLcrM in (l-l-l)-dimensions; Starting from the Haldane 
conjecture, it is by now well-known that the different 
nature of the ground states of the S = integer and 
S ~ half-odd integer spin chain is reflected to the value 
of the 6'-term, 6* = 27r x 5 (mod 27r) in the 0(3) NLcrM. 
Quite surprisingly, in (l+l)-dimensions, the integrabil- 
ity of the 0(3) NLaM allows us to "prove" that the 0(3) 
NLcrM &t 9 ~ TT flows, in the RG sense, into a gapless 
critical point (the SU(2) Wess-Zumino-Witten theory at 
level one). The low-lying excitation at the critical point 
consists, not of magnons, but (fermionic) spinous. When 
detuned from 9 = n (by breaking a link parity symmetry. 
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say, with bond dimerization in the Heisenberg exchange 
couphng), the ^-term flows to or 27r, where a system is 
in a gapped paramagnetic phase. 

We now present some speculations on the effect of the 
topological term in (2+1) dimensions in the region when 
the fluctuations of the NLcrM order parameter is strong. 
Let us start from the ordered phase of the NLcrM; it ex- 
ists for the small NLcrM coupling, i ^ 1 , and this should 
be true for any value of 6*. As we crank up the NLcrM 
coupling constant t, the interactions between "magnons" 
become stronger. It would then seem reasonable to as- 
sume that, even in the presence of non- vanishing 9, we 
reach a critical point by eventually destroying the (long- 
range) order, at which the NLcrM order parameter de- 
velops a critical (power-law) correlation with vanishing 
expectation value, v = 0. 

Recall that in the conventional NLcrM without any 
topological term, the nature of such transition can be 
captured by e.g. the 2 + e expansion. On the other 
hand if we stay exactly in three space-time dimensions, 
which is necessary if we are after the effect of the topo- 
logical term, the NLcrM is not perturbatively renormal- 
izable (although it may be renormalizable in terms of 
the non-perturbative RG). It should the be considered 
as an effective field theory, and is not UV complete. As 
we approach the putative quantum critical region (the 
strong coupling region of the NLcrM) described above, 
the NLcrM description should be replaced by some other 
description. The question then is, to which extent we can 
deduce such UV (high-energy) description for the case of 
non- vanishing 9. 

In terms of the corresponding fermionic model Sec. 



II B , which should be regarded as a more microscopic 
(or UV) description, increasing t (the coupling constant 
of the NLcrM) corresponds to reducing the four fermion 
coupling constant g. In terms of the fermionic phase di- 
agram (Fig. |l|), we thus approach the phase boundary 
from the ordered side. In particular, \v\ — corresponds 
to the case of = 0, and in this case, we approach the GN 
critical point separating the ordered phase and paramag- 
netic phase (Dirac scmi-mctal phase). One would then 
conclude that a UV critical point of the NLcrM at 9 = tt, 
if exists, would correspond to the GN fixed point. 

This is a tempting, but of course very dangerous, argu- 
ment. We would, however, expect, by matching with the 
fermionic description, the following: at the UV critical 
point of the NLcrM at 9 = tt, not only the NLcrM order 
parameter develops critical correlation, but its dual or- 
der parameter also does so. (These order parameters are 
called vi and V2 in the beginning of this section). This is. 



of course, a complete surprise from the NLcrM point of 
view: There is no inkling of V2 order parameter what so 
ever in the ordered phase, and controlling coupling con- 
stant of the NLcrM {t in Eq. ^3|) has nothing to do with V2 
order. Nevertheless because of the 6'-tcrm, making t large 
(destabilizing the order parameter vi) somehow makes V2 
correlation stronger, since the V2 order parameter within 
the NLcrM is realized as a defect in the vi parameter. 
This is in fact in line with the non-trivial identification 
of the VBS order parameter as the skyrmion creation op- 
erator in the (2-1-1) dimensiavial 0(3) NLcrM augmented 
with the Berry phase terniBilE] The VBS phase arises as 
a proliferation of monopoles in the presence of non-trivial 
Berry phases. Thus, while the more precise nature of the 
putative quantum critical point at = tt is difficult to 
study, overall physical picture related to the physics of 
0-term seems to be inferred from the fermionic GN model 
description. 

Another interesting issue is the RG flow which incor- 
porates both t and 9 (the couplings of kinetic and topo- 
logical terms in the NLcrM). In the weakly coupled re- 
gion {t/A ^ 1), the presence of the non- vanishing 9 little 
affects the RG flow, while this may not be the case in 
the strongly coupled region {t/A ~ 1). In fact, this is the 
case for the NLcrM on the Grassmannian manifold in two 
dimensions, which is relevant to the-jtwo parameter scal- 
ing flow of the quantum Hall effect .cil Again by matching 
the NLcrM description with the fermionic GN model de- 
scription, since the 0-term is correlated to the mass term 
mtj-tpTit^' (sec Eq. |4^), it would be then interesting to 
"guess" the RG flow of the 6'-term from the RG flow of 
the mass term. Within the large-Nf or e expansion, the 
scaling dimension of the fermion mass rubip'^bil' the GN 
critical point is very close to the scaling dimension at the 
trivial (non-interacting Dirac) fixed point. Thus, when 
l/Nf or e is small enough, the mass term is a relevant 
perturbation to the GN fixed point, and it grows under 
the RG. This would suggest the 6'-tcrm is also relevant 
in the strong coupling region of the NLcrM at 9 = tt; the 
deviation of the 0-angle from 9 = tt grows under the RG. 
However, in principle, more exotic possibility, i.e., the 9- 
term at the non-trivial critical point being irrelevant, can 
also be realized. We leave these issues for future studies. 
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